A New Limit on the Distances of Nuclei UHECRs Sources 



O 
(N 



00 



Tsvi Piran^'B 
' Racah Institute of Physics, The Hebrew University, Jerusalem, Israel 

Recent evidence from the Pierre Auger Observatory suggests a transition, at 5 EeV-lOEeV in the 
composition of Ultra High Energy Cosmic Rays (UHECRs), from protons to heavier nuclei such as 
iron. I consider here the implications of the heavier composition on the sources of UHECRs. The 
smaller magnetic rigidity implies that nuclei UHECRs are: (i) More easily accelerated (ii) Local, as 
they can diffuse only a few Mpc from their sources before disintegrating (iii) Isotropic, because large 
deflections in the extra galactic and the galactic magnetic flelds erase the directional information. 
Uncertainties in the strength and structure of the the extra galactic magnetic field (EGMF) makes it 
difficult to estimate the overall effects. However, with typical reasonable parameters of a few nC and 
a coherence distance of a Mpc the distance a nuclei UHECR above the GZK energy traverses before 
photodisintegrating is only a few Mpc. In spite of the significantly weaker limits on the luminosity, 
Cen A is the only currently active potential source of nuclei UHECRs within this distance. The large 
deflections erases the directional anisotropy expected from a single source. If indeed the composition 
of above-GZK-UHECRs is iron and if the EGMF is not too small then Cen A is the dominant source 
of observed nuclei UHECRs above the GZK limit. 

PACS numbers: 96.50.S-, 13.85.Tp, 98.70.Sa,98,62En,98,70Sa 
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I. INTRODUCTION 

Recent observations of the Pierre Auger Observatory 
(PAO) show a transition in (Xmax) and in _RMS'(X„iax) 
between 5EeV and lOEeV [l[. These transitions are in- 
terpreted as reflecting a transition in the composition of 
Ultra High energy Cosmic Rays (UHECR) in this energy 
range from protons (denoted hereafter as pUHECRs) to 
intermediate weight nuclei and in particular towards iron 
(denoted hereafter nUHECRs). Like protons, iron nu- 
clei also suffer a strong annihilation above the Greisen- 
Zatsepin-Kuzmin (GZK) 0, i] energy of ~ 50EeV. How- 
ever, the larger charge and mass affect both the accel- 
eration and the propagation of nUHECRs leading to a 
drastically different view on their possible sources. 

I begin by a brief review of the observations, (i) Both 
PAO [4] and HiRes [51 found a decline in the spectral slop 
of UHECRs above ~ 50EeV. This is roughly at the en- 
ergy for which the Greisen-Zatsepin-Kuzmin (GZK) [3,y| 
suppresion due to the interaction of the UHECRs with 
the CMB is expected. The relevant interactions are pho- 
topion production for protons and photodisintegration 
for nuclei, however coincidentally, the threshold for both 
interaction is roughly at the same energy, (ii) As al- 
ready mentioned, the PAO reports a transition in com- 
position from protons to intermediate mass nuclei, more 
specifically iron, between 5EeV and lOEeV [1]. These 
results are consistent with a recent analysis of Yakustk 
data [6| but are not supported by HiRes that finds a pro- 
tonic composition (7|. (iii) The overall sky distribution 
of the arrival directions of UHECRs below the GZK en- 
ergy is isotropic (see, e.g., |8l-[l0|). However, at higher 
energies, which are of more interest to us, analyses of 



different experiments lead to somewhat conflicting con- 
clusions. For example, AGASA finds an excessive num- 
ber of pairs and one triplet in the arrival direction of CRs 
above 40EeV suggesting small scale clustering [ll|. On 
the other hand, the HiRes stereo data are consistent with 
the hypothesis of null clustering J^j while the autocor- 
relation analysis of the PAO data reported a weak excess 
of pairs for E > 57EeV [ll]. The PAO data shows a 
large scale isotropy [iJl however, PAO [lj| found a cor- 
relation (within a radius of 3.1°) between events above 
57EeV with AGNs located closer than 75 Mpc in the 
Veron-Cetty & Veron catalog [l^- The HiRes data, 
however, do not show a correlation of the highest en- 
ergy events with nearby AGNs [16| , but there have been 
controvercial claims of correlation of the HiRes data with 
distant BL Lac objects [17|. Finally, the PAO also finds 
for E > 57EeV a "hot spot" in the direction of Cen A 
[18| . Different authors make different cuts and attribute 
different number of UHECRs and a different statistical 
significance to this "hot spot". 

I consider here the implications of an iron composition 
of the highest energy CRs on the acceleration (S|TT| and 
on the propagation ( mil^ of UHECRs. I examine the 
implications to relevant sources in mV\ 



II. ACCELERATION 

The smaller magnetic rigidity of the iron leads to 
several important differences on the nature of possible 
sources of nUHECRs as compared with sources of pUHE- 
CRs. First, it eases the strict acceleration constraint [IQi]: 



RB > 10^^ (cm G) E2oT^(3-^Z. 
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where T and /3 are the Lorentz factor and velocity of the 
source. Generally Q^ = Q/W and E2Q = E/10'^°eV but 



for the electric charge ^26 = ^/26, in units of the elec- 
tron's charge. This criterion sets a limit on the Poynting 
flux luminosity of the source (see e.g. [20|)' 

L > 1.5 X 10^2 (erg/s) {T^ / fi)El^Z^;^ . (2) 

A comparison of synchrotron losses with the acceleration 
rate limits the source's magnetic field: 



B < 3 G r^S^O (^26M56)"^^2 



(3) 



where A^^, is the atomic weight in units of the proton's 
mass. An additional condition on the magnetic field 
arises for nuclei. The energy of a synchrotron photon 
generated by the gyro motion of the nucleolus shouldn't 
be large enough to destroy the nucleolus: 



B < 0.3 G Z2QA-^{hv synch /^Me\). 



(4) 



Both limits on the magnetic field in the source (Eqs. [3] 
andHI are easily satisfied for AGNs. A related condition 
arises on the strength of the radiation field at photon 
energies of a few MeV within the acceleration regions. 

As can be intuitively realized the acceleration of a par- 
ticle with a charge Z is much easier than the acceleration 
of a proton. These eases the conditions at the accel- 
eration site. Particularly important is the much lower 
(ex Z^^) limit on the sources' synchrotron luminosity. 
Eq. [2] implies a well known drastic difference between the 
possible sources of nUHECRs and of pUHECRs. While 
AGNs with L > 10'*^erg/s are rare and none exists 
at present within the GZK distance Lqzk ~ lOOMpc, 
AGNs with L > iC'^erg/s are numerous. The lower 
rigidity relaxes the most critical GZK problem, the lack 
of suitable accelerators within the GZK distance. 



III. PROPAGATION 

Not less drastic effects arise concerning the propaga- 
tion of the nuclei from the source to Earth. The lower 
rigidity leads to a much smaller Larmor radius: 



i?L = 4 Mpc E2oB-^Z^^^ , 



(5) 



where B is EGMF. The expected For Rl < A, the maxi- 
mal coherence length of the EGMF, the particles diffuse 
in Kolmogorov regime with a mean free pass, I: 

I « (i?LA2)i/^ = 1.6 Mpc {E2o/BnGZ26)^^^X%lc: (6) 

where the expected value of A, is of order 0.1 — lOMpc. 
The corresponding diffusion coefficient D{E, B) satisfies: 

D « 0.85 MpcVMy (£;2o/^26S„g)'/^A^/^^. (7) 

The most energetic nUHECRs may have Rl{E)>\ and 
could be in the transition region between the Kolmogorov 
and the Bohm regimes. For simplicity I assume in the 



following that all relevant particles are in the Kolmogorov 
regime. 

The maximal distance, dmax{E), that an average 
nUHECR with energy E traverses before photodisinte- 
grating is: 



^m.n.T, — 
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Mpc 



IGGMpc \ B„G^26 



(8) 

With typical parameters dmax ^ Lqzk- This implies 
that nUHECR sources are much nearer than what was 
expected earlier. Eq. [8] is the second drastically different 
feature of a nUHECR as compared with a pUHECR. 

Consider now a source at a distance d from Earth 
that emits nUHECRs from time Ton until Toff. The 
maximal propagation time of the CRs is tmax = 
Min(L(32/f /c, Ton). There is a negligible contribu- 
tion from times prior to the arrival of the diffusion 
front [2l|, so the minimal propagation time is tmin = 
Ma,x{d'^ /3D,Toff). For a uniform emission with a con- 
stant rate no the particles distribution satisfies: 



n{r,t) 



2no(t, 



-1/2 
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^■kD{E,B)Y/^ 



(9) 



If the original spectrum of the particles at the source 
is E'P the observed spectrum will be E'^'^p^^^^'^) where 
the diffusion coefficient satisfies D ex E^ . A consistent 
solution with the observed E~'^-'^ spectrum is obtained 
for an injection spectrum of E^'^-'^ if the nUHECRs are 
in the Kolmogorov diffusion regime for which /3 = 1/3. 
The spectrum might bc^ too steep in the Bohm regime for 
which /3 — 1. Kolmogorov diffusion holds if BnG^Mpc<^5. 
Given the uncertainties in the EGMF structure this is a 
viable possibility. 

The CR distribution will be anisotropic with [21|: 



/(6',r,i) = (1 + a cos 61) 



n{r,t)c 
An 



(10) 



The anisotropy in the CR distribution is the flux 
weighted average of d/2tc so a « d/{&tminc)- For a 
source that is still active now: 



«.^.0.25 
2dc 



E20^i\ 



1/3 , 

Mpc \ / 3.8Mpc 



BnGZ26 



(11) 



For a single source, given the observed isotropy of the 
UHECR [ii, [Hi Eq. [TT] might seem to be a problem. 
However, for lower energy UHECRs below the GZK en- 
ergy the diffusion coefficient is much smaller and the 
anisotropy due to a local source will be erase. Above 
the GZk energy this anisotropy will be smaller if the 



^ The observed spectrum above the GZK energy is steeper than 
-2.7, but this is attributed to the GZK suppression. 



source has turned off at Toff > (P/3D. Additionally, 
and more important, this (energy dependent) anisotropy 
will be erased by strong deflections in the Galactic mag- 
netic fields. The magnetic field in the Galaxy is of or- 
der 4/iG. It is composed of an ordered component and 
a random component. The Larmor radius in this field 
of even the highest energy nUHECR is i?Lg^,<kpc. For 
most arrival directions this is smaller than the distance 
a nUHECRs traverses in the Galaxy leading to a signifi- 
cant defiections that will erase the anisotropy of the fiux 
reaching the Galaxy. Numerical simulations of pUHE- 
CRs find Galactic defiections angels of a few degrees [23] 
corresponding to large (90°) defiections for nUHECRs. 



IV. SOURCES 

The new limits on nUHECR sources are quite differ- 
ent from those for pUHECRs. (i) The lower luminosity 
limit (Eq. [2]) increases significantly the number of source 
candidates, (ii) On the other hand the small maximal 
distance (Eq. El) decreases significantly the volume where 
the sources can be. (iii) The larger deflections in both 
the EGMF and in the Galactic magnetic field erase most 
if not all directional information on the source. 

Clearly if the EGMF is suflficiently small B„g < 0.01, 
as has been recently suggested [2J, l2a| , its effect on the 
nUHECRs would be minimal and there won't be a signifi- 
cant new lower limit on the distance of nUHECR sources. 
However, (i) and (iii) would still be applicable. The later 
will be valid since nUHECRs will still suffer strong de- 
flections in the Galactic magnetic held. Thus even in 
this case a significant correlation of the arrival direction 
of nUHECRs with their sources should not be expected. 

However, the common understanding is that the 
EGMF is stronger, of the order of a few nG. Upper limits 
on rotation measures (RM) of radio signals from distance 
Quasars were used [23 to set a limit of B < lO^^^G for 
a homogenous field and a B <nG for a random field with 
a coherence scale Xmpc ~ 1- However Qb — 1 was used 
in these estimates to obtain the electrons' density. With 
Hb ~ 0.04 the actual upper limits are larger by a fac- 
tor of 20, reaching Bnc ~ 20 [2gl. More recent work 
suggest that the EGMF follows the large scale structure 
with magnetic fields of order a fiG within clusters |29| 
and 10 — lOOnG within superclusters [30|, |3l|. There 
has been an extensive work on propagation of UHECRs 
in such magnetic fields (see e.g. [S^, [33] and references 
therein). However, most if not all, this work is concerned 
with the propagation of pUHECRs. Scaling the results 
to nuclei suggests that nUHECRs are in the strong scat- 
tering regime. For example [32] find, for a realistic mag- 
netic field distribution that follows the large scale struc- 
ture, that the optical depth for a magnetic deflection of a 
pUHECR is larger than unity and the overall deflection 
angles are of a few degrees. These results suggest that 
the local universe is magnetically opaque for nUHECRs 
whose typical deflections are of order unity. 



If the EGMF is of order of a few nG or if the magnetic 
field in the Virgo supercluster is of order O.l/iG then 
it follows from Eq. [8] that the dominant source is the 
nearest active radio galaxy Centaurus A (Cen A). M87 
may possibly add a minor contribution. At a distance of 
only 3.8 ± O.lMpc [3J] Cen A proves to be an excellent 
source candidate. Photon emission from the nucleus of 
the galaxy has been detected in the radio, infra-red (IR), 
X-ray, and in the GeV-TeV range. Radio and GeV emis- 
sion was also observed from the large radio lobes that 
extend up to 250 kpc. Cen A is also the nearest extra 
galactic source of TeV photons which, apart from UHE- 
CRs, are the highest energy particles observed so far on 
Earth. As such they may be a good indication for pro- 
duction of UHECRs as well. In the past, Cen A was 
already suggested as a major (possibly only) source of 
UHECRs 0,|3i-[33. This possibility has received a lot 
of attention recently following the observations of a pos- 
sible (^ 2% significance) concentration of UHECR with 
arrival directions within a few degrees from Cen A |18| . 
However, if UHECRs are indeed iron nuclei this associa- 
tion is most likely spurious (unless both the EGMF and 
the Galactic magnetic fields in the direction of Cen A 
are particularly low). The role of Cen A as the dominant 
nUHECR source arises from the small maximal distance 
that nUHECR can propagate even in a modest EGMF! 
Cen A is the strongest and possibly only known source 
within the reasonable distance. M87, at 16Mpc, might 
also contribute, if the EGMF is on the lower side. 

The observed present total luminosity of Cen A is 
^ 10'''^erg/s, of which about half is in high energy [38{: 
comfortably above the synchrotron luminosity limit of 
Eq. [5] Let ejj be the efficiency of UHECR production 
compared to photon production extrapolated to lOEeV 
using equal power per decade and e^ be the fraction of 
energy that goes to accelerate iron nuclei. If Cen A is the 
only source in the magnetic GZK volume of {4:Tr/3)d^^^ 
then the local UHECR injection rate is £uei5.5 x 10^^ 
erg/Mpc'^/yr. This is easily consistent with the observed 
value of 0.45 ±0.15 x 10''^ erg/Mpc^/yr [H, even consid- 
ering the fact that the majority of the accelerated par- 
ticles are lower energy protons and only a small fraction 
of the total UHECR energy might go to iron nuclei. 

Cen A is a nearby source (in terms of d/l). This makes 
the analysis somewhat problematic in view of the large 
uncertainty in the value of the intervening magnetic field 
and hence in the mean free path, I, and the correspond- 
ing diffusion coefficient. With a low value of the magnetic 
field the highest energy nUHECRs could be in the Bohm 
diffusion regime and with d ^ I. This would, of course 
make some of our simple estimates, in particular Eqs. |9l 
[TUl and [TT] for the nUHECR density and the anisotropy 
invalid. Another complication, which is ignored here, 
arises if there are large voids (of orders of tens of Mpcs) 
in the magnetic fleld structure that allow nUHECRs to 
propagate in straight lines over large distances. Detailed 
simulation 40;] suggest isotropy for pUHECRs emitted 
from Cen A with B ~ 1/iG. This corresponds to 50nG 



for nUHECRs. This is comfortably within the expected 
estimates within the local supercluster. Similarly, sim- 
ulations of pUHECR propagation in the Galactic fields 
|23| yield deflections of a few degrees, which imply de- 
fections of order unity for nUHECRs. It is interesting to 
note that within this model one cannot expect a com- 
plete isotropy as both the EGMF magnetic fields and 
the galactic field could induce an anisotropy arising from 
their own structure. For example a strong EGMF in the 
Virgo supercluster could induce excess of particles around 
the super-galactic plane. 



V. CONCLUSION 

// UHECRs above the GZK are dominantly iron nu- 
clei and if the EGMF is larger than a Jew nC then the 



largest possible distance of the source would be only a 
few Mpc. Even though it is much easier to accelerate 
iron nuclei than protons the only available source within 
such a distance is Cen A. Large deflections of order unity 
in both the EGMF and in the Galactic magnetic fields 
would erase, in such a case, most the directional infor- 
mation on the UHECR sources. 
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